Intermediate spin approach to the tunneling of nano-magnets. 
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The theory for the quantum tunneling of nano-magnets is developed within the intermediate spin 
framework. Periodic magnetic effects are seen to reflect that associated with a change of flux by a 
single flux quantum $o- Essential are Schrodinger cat wave functions which involve superpositions 
of different magnitudes of the applied magnetic field. For systems in which the tunneling paths 
are not co-planar the theory leads to essentially complex magnetic fields although the expectation 
values of the fields remain real. In general the ground states correspond to minimal uncertainty 
squeezed states. The degree of squeezing depends on the anisotropy parameters. 
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I. INTRODUCTION 

The quantum properties of very small magnets often 
called nano-magnets have generated a good deal of in- 
terest over the past decade. Much of the activity has 
focused upon the demonstration of quantum mechanical 
tunneling in its various aspects |J. Here attention is 
drawn to another purely quantum phenomenon namely 
intermediate spin [§,[||. Intermediate spin comes as part 
of the baggage of intermediate statistics [||, i.e., anyons 
and two dimensions. The intermediate statistics and spin 
reflect a statistical parameter a s which is conventionally 
a s = for bosons and whole integer spin and a s = 1 
for fermions and half-integer spin. It does not seem to 
be generally appreciated that there are a good number 
of (real and numerical) experimental manifestations of 
intermediate spin in relatively simple magnet systems in 
a magnetic field. The principal purpose of this paper 
is develop a formalism for quantum tunneling of nano- 
magnets based upon the theory of intermediate spin. 

It will be show that a magnetic field directed along a 
suitable symmetry axis can cause, e.g., a physical whole 
integer nano-magnet to have an low energy level structure 
which is identical to the equivalent half-integer system, 
this including a Kramer's degeneracy. Zero and a certain 




FIG. 1. A particle of charge q is confined to a unit circle. 
A tunneling problem might be defined by adding a potential 
V such that V{<j>) = V^ + n). If this potential has minima at 
= and 7i and barrier between is higher than the zero point 
energy then there will be a ground doublet with a splitting 
A. A perpendicular field exerts no force on the particle but 
does change the slitting as explained in the text. 



fields lying between the half-integer values correspond to 
whole integer points. Intermediate fields correspond to 
intermediate spin. 

The suppression of the ground doublet tunnel split- 
ting which necessarily occurs at half-integer points has 
been previously recognized by Garg Q. What has been 
labeled "topological quenching" is in fact just a pseudo- 
Kramer's degeneracy of the ground state. The implica- 
tions of intermediate spin are much stronger since all low 
lying levels must be degenerate at exactly the same field. 

It is a commonplace observation, for the situation illus- 
trated in Fig. (1), that magnetic effects are periodic with 
period $0, the flux quantum. For the problem explained 
in the figure caption in particular the tunnel splitting, 



A = A |cos27T$/$o|- 



(1) 



For the magnetic systems the symmetry which will be 
considered in most detail here is that appropriate p,p| to 
Feg, Fig. (2a). This corresponds to the most general 




FIG. 2. The length of the axis indicates the energy gain 
when the spin S lies along that direction, thus (a) for the sym- 
metry appropriate to Fes when all three axes are inequivalent 
the vertical ?/-axis is easy while the field is directed along the 
hardest x-direction. For (b) appropriate to Mni 2 with equiva- 
lent x and y directions, the z-direction is again easy. The field 
is along the hardest direction in the x — ^/-plane. For cubic 
symmetry, illustrated is the possibility with the field along a 
hard four fold axis. The perpendicular axes directed to the 
corners of the cube are all equivalent easy directions (in the 
absence of a field) . Intermediate spin effects also occur if the 
field is directed along a three fold axis. 



quadratic nano-magnet Hamiltonian: 

H = -DS Z 2 + E [S x 2 - S y 2 ] 



(2) 



which without loss of generality has two positive 
anisotropy parameters E < D, see the Appendix. A mag- 
netic field is directed along the hardest, i.e., x-direction. 
This is fully equivalent to the particle problem of Fig. (1) 
with $/$o = gii B H/Ah 2 ^2E{D + E). The agreement 
with Eqn. ([!]) is illustrated by Fig. (3). The deviations 
from the | cos | behavior for larger fields is due to the fi- 
nite value of S. It is larger for the value S — 10 relevant 
for Feg. 

Under an adiabatic change of magnetic field by $0 the 
S z quantum number m — > m + 1 for the angular momen- 
tum quantized along the magnetic field direction. What- 
ever the mixture of states with different m, the period 
measured in terms of the magnetization along the field 
direction is automatically %/i#. However it is a trivial 
consequence of time reversal invariance that the states 
with fluxes ±<I> are degenerate and so there is the possibil- 
ity to construct "Schrodinger cat" states which comprise 
a linear combination of states with the macroscopic field 
H in both its positive and negative directions. In this 
case the magnetization period need no longer be %/i#. 
It will be shown, with this symmetry, such a combina- 
tion exists in all energy eigenstates and that the period 
is indeed a fraction of hgfiB • 

An even more bizarre phenomenon occur when either 
the field is not directed exactly along the appropriate 
symmetry axis and/or if the least action tunneling paths 
are not co-planar. For these cases the energy eigen- 
states correspond to suitably weighted "Schrodinger cat" 
combinations with oppositely directed complex magnetic 
fields. The complex weighting is such that the expec- 
tation value of the physical magnetic field remains real. 
Thus, again for the Fes symmetry, a field H t which is 
perpendicular to both the hard and easy directions cor- 
responds to a pure imaginary field iH t directed along 
the symmetry axis, i.e., the tunnel splitting is given 
by Eqn. @) with $/$ = igfi B H t /4,h 2 y/2E(D + E) 
so that cos — > cosh. In the presence of both fields 
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FIG. 4. To the right the thick line is approximation 
Eqn. (1) with a complex field B c — B(cos6 + isinO). The 
thinner line corresponds to exact diagonalization. To the left 
is experiment pi with a logarithmic scale. The imaginary part 
of the flux causes the tunnel splitting to increase with field. 



$/$o = gfJL B (H + iH t )/4h 2 v /2E(D + E). The result- 
ing approximation is compared with the exact result and 
experiment [|6| in Fig. (4). 

For the much studied |7|,|8[] molecular magnet Mni2 it is 
impossible to not have a complex field. With the physical 
field in the direction indicated in Fig. (2b) the resulting 
complex field has equal real and imaginary parts in the 
large S limit, see Fig. (5). There are two equivalent tun- 
neling planes which both make an angle = 7r/4 with 
the applied field (and 7r/2 with each other). The wave 
function for the planes is added so that there are two sets 
of paths, one from each plane, which are equivalent and 
which interfere. The resulting expression for the tunnel 
splitting is shown in the figure caption. As the figure 
illustrates this splitting increases very rapidly with field 
and by several orders for the S oscillations. 

The last important principle is illustrated by applying 
a field in the hardest direction for a system with cubic 
symmetry, Fig. (lc). This is a four fold axis and group 
theory insists that the lowest tunnel split multiplet com- 
prises four states. For whole integer S it turns out that 
there are two singlets and one doublet, however for a half- 
integer point this must reduce to two doublets, reflecting 
the pseudo-Kramer's theorem. That this is the case for 
large S is illustrated by the exact results of Fig. (6). 

Classically for magnet such as Fes the equilibrium 



Splitting 




FIG. 3. Properties are periodic in the field h = H/g/^sh for 
S — 35. The deviations from A = Ao| cos27r$/$o| become 
smaller as S increases. 



, Field 

2 4 6 S 

FIG. 5. The complex field result is 

A | cosh(7r/i/p) cos(tt/i/p)| where p « 2S 1/2 (12aD 3 ) 1/4 . This 
corresponds to a complex magnetic field in which the real and 
imaginary parts are equal. This exact result is indicated. The 
difference becomes smaller with increasing S. 
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FIG. 6. For the four fold cubic axis there are four tunnel 
split levels. Shown are the splittings between levels 1 and 
2 and between 3 and 4. The latter have been moved up by 
0.001. Illustrated is the fact that both splittings are zero 
for the same field. This reflects a pseudo-Karmer's theorem 
for half-integer points. (Notice that the splitting increases 
corresponding to a complex magnetic field.) 

magnetization lies along its easy axis and is maximal. 
The existence of quantum fluctuations is reflected in the 
fact that the expectation value of this easy axis magne- 
tization is hS and not h^S(S + 1). In the absence of 
a hard direction in the plane perpendicular to the easy 
axis, the quantum fluctuations measured by (S x ) and 
(S x ) ~ h 2 S/\/2 are equal and in fact minimal. In 
the terminology used for fluctuations in quantum op- 
tics, these are minimal un-squeezed states. Adding an 
anisotropy energy E, which creates a hard axis in this 
plane as is appropriate to Feg, squeezes the quantum 
fluctuations so that (S x 2 ) ^ (S y 2 ) but leaves minimal 
states. The present approach thereby lends itself to a 
natural definition of squeezed spin states || . 

Apart from the simple periodicity [pi of Eqn. ffl) the 
intermediate spin effects described above have yet to be 
experimentally verified. 



II. INTERMEDIATE SPIN 

Returning to the basic charged particle problem of 
Fig. (1), the relevant group is SO (2) which is equiva- 
lent to U(l), i.e., is abelian with a single operator S z , 
the generator of rotations e %Sz $ \ where the z-axis is per- 
pendicular to the circle and where <\> is the angle to the 
x-axis. For a unit circle S z = p^ the momentum conju- 
gate to (/>, i.e., [4>,Pcf)] = ih. The eigenstates of S z and 



p^ are VVnO) 



but for particles on a cir- 



(2tt) 1 /2 c 

cle there is no reason it insist on m being an integer 
reflecting the fact that the spectrum of S z for SO (2) is 
continuous. In general for this singular gauge i/j m ((f)) is 
mult i- valued. A convenient single valued function can be 
defined by taking some convention for ip m (<fi) in the inter- 
val {0, 2tt} and then analytically continuing this function 
to the interval { — oo,+oo}. For any physical potential 
V{<j>) = V(<f> + 27m); n = 0, ±1, ±2 . . ., and implies that 
only the basis states ^^s +n (0); n — 0, ±1, ±2 . . ., mix. 



It follows for a given statistical parameter a s the energy 
spectrum is discrete. This can be stated differently: Since 
V((f)) is periodic the solutions obey Floquet's (Bloch's) 
theorem, i.e., ipk(<ft) = e %k ^ Uk{(j)) where Uk(<fi) is periodic. 
The reciprocal lattice vector K = 2tt /(a = 2tt) = 1 and 
the energy e(fc) = e(fc + nK); n = 0,±1,±2... i.e., 
is periodic. An examination of the Fourier components 
immediately implies k = a s /2. Given that V((j)) has a 
point with reflection symmetry the ip±k(<ft) are degener- 
ate and the most general solution for a given energy is 
Aipkify+Bip-kifi) for which the spectrum of 771 = n±^f , 
i.e., there are twice the number of observable S z values 
except importantly for the whole and half-integer points 
when a s = or 1. This more general method of extrapo- 
lating between whole and half-integer spin which will be 
found applicable here. 

In the non-singular gauge the basis set ip n (^>) nas 
n — ±1, ±2 . . . and the charged particle encircles a mag- 
netic flux <£ = a s $>o/2 where <l>o is the appropriate flux 
quantum. For the same V{4>) there exists an eigenstate 
Uk((f>) with the same energy as in the singular gauge, 
i.e., the energy spectrum are identical and two prob- 
lems can mapped to each other and both reflect inter- 
mediate spin. In the non-singular gauge the combina- 
tion Auk(4>) + Bu-k(4>) comprises a linear combination 
of states with the flux <£ in both its positive and negative 
senses. Evidently the expectation value of the flux (3>) 
is in general smaller than a s <&o/2. This possibility of a 
linear combination of different flux states does not seem 
to have been envisaged in connection with intermediate 
statistics although it would arise if, e.g., quantum coher- 
ence in SQUIDS could be observed. It is realized in the 
systems, e.g., Feg, discussed here. 

III. SPIN MODELS 

It has become usual to re-write Eqn. (2) with the x 
and z-axes interchanged. Thus with fields the general 
quadratic nano-magnet Hamiltonian becomes, 



H = -(D-E)S X 2 + 2ES Z 2 
-ES 2 - hS z 



hiS x — hfSy, 



(3) 



see the Appendix. Unless stated otherwise it will be 
assumed that the physical spin value S is whole inte- 
ger. The value S is that which would be assigned by 
the usual arguments, e.g., for a set of ferromagnetically 
ordered spins this is simply the sum of the spin values. 
The field h is parallel to the hardest and transverse to 
the easiest axis. It is equivalent to that introduced in 
the non-singular gauge. The other transverse h t and the 
longitudinal hi fields simply modify the potential V. 
Since the length of S is fixed, classically the energy, 

E(0, (j)) = h 2 S 2 sin 2 6 [-(£> - E) cos 2 <\> + 2£sin 2 0] 
- hhS cos - hihSsmOcoscj)- h t hS sin sin (j) - Eh 2 S 2 , (4) 



is a function of the usual two angles and 0, see Fig. (7). 
Quantum mechanically Eqn. (|4|) is the energy expec- 
tation value for the coherent state \Q,(/)) obtained by 
rotating \S Z = S) through 9, (j). It is these coherent 
states which appear in the usual functional integral for- 
mulation |l|. However the \6,<j)) are very much over- 
complete. In fact the states |</>) = \6 = §,</>) = (25 + 

D 1/2 EL-5 



1( ^|m), where \m) = \S Z = m), 



remain 



over-complete. If 



E 



< a(m)\m) then ip(<f>) 



m) = e; 



e* m ^a(m), i.e., is a Fourier transform of 



the m-space wave function a(m). The function ^((j)) is 
defined on the interval {— oo, +00} and the inverse trans- 
form ip(m) = ± /+~ d</nl>(<l>) = Y,m> S ( m ~ rn')a(m') 
so that \ip) = J_°° dmifj(m)\m), i.e., i/j(m) is a defini- 
tion of the m-space wave function appropriate when m 
is considered as a continuous variable. Clearly S z ip((j)) = 
J-™ dme l7nc ^hmij)(m) = p^ip^cf)) where p^ = —ihjz is 
the same as that defined earlier. Also evident is that 
multiplication by e %< ^ amounts to a displacement of m by 
unity and so, e.g., S+ = e^[h 2 S{S + 1) -ityftty + 1)] 1/2 . 



IV. TOPOLOGICAL CONSIDERATIONS 

In terms of the unit sphere, the system of coordinates 
(j),P(f) might be considered as a particular projection in 
which the axis of quantization plays the role of the North- 
South axis as shown in Fig. (7). In these terms, the angle 
(j) is evidently the longitude while p^ is the projection on 
the z-axis and is therefore the sine of the latitude. That 
these are canonical coordinates and momentum empha- 
sizes the evident fact that the unit sphere is the phase 
rather than simply the coordinate space. 



In order to calculate the Berry phase |10| a unit vector 
R is directed along the easy, i.e., x-axis for Eqn. (||) (the 
z-axis for Eqn. ([!])). It is particularly easy to determine 
this phase for an adiabatic process in which the system 
is rotated about the z-axis (the x-axis for Eqn. (El)). For 
the singular gauge, if the wave function comprises a single 
Floquet (Bloch) wave e %k ^ Uk{(j)) this is trivially 



G B = 2tt/c = 2tt— 
2 



(5) 



which is the obvious generalization of the result QS where 
17, here 27r, is the solid angle subtended by the path of 
R. However for the combination Ai/jk(<f>) + Bifj-k((j)) the 
result is tanB B = [(A - B)/(A + £)]tan27r^f which 
agrees with Q 9 ^- only at integer and half-integer points. 



V. PSEUDO KRAMER'S OPERATORS 

If it is possible to transmute a whole into a half-integer 
spin then at the half-integer points there should exist an 
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FIG. 7. Mapping of the unit sphere to a two dimensional 
rectangle 

operator K, the equivalent of the time reversal operator 
for h = 0. This must commute with Ti and if \E n ) is an 
eigenstate then K\E n ) should be orthogonal (or zero). 
For the current model Preda and Barnes ]ll| have con- 
structed such operators. 

The simplest case is with hg = ht = and h = 
h^/2E(D + E). The operator K t = Tr where T is the 
conventional time reversal operator and where 



r = VD + ES X + iV2ES y . 

It is easy to show that with this field the negative semi- 
definite, 



n 



-(Tr) 



so that trivially [H, K\\ = independent of S and the 
parameters D and E. As will be shown in Sec. VII, 
with hi = h t = there is an eigenstate of the form 
\ E n) = S p a ( 2 P)l 2 P)> where S z \m) = ftra|ra), i.e., only 
states with even values of m are present. Since r only 
has matrix elements between even and odd such states it 
follows that K\\E n ) is orthogonal to \E n ) and hence for 
this value of the field all states are double degenerate or 
K\\E n ) = 0. Clearly this latter state has energy E n = 0. 
It is evident that when E = there is only one such 
state and by a continuity argument this must be the case 
for all E. It follows that there are S degenerate pairs 
when h = y/2E(D + E) fully mimicking a half-integer 
spin except for the highest energy singlet. 

It is therefore a remarkable fact that this degeneracy 
for fields which are exact multiples of h^2E(D + E) ex- 
ists for all but the highest energy states for arbitrary 
values of integer S and independent of the existence of 
a tunneling barrier, etc. For example it is possible to 
explicitly construct the S — 1 matrices and show that 



for ft = h^/2E(D + E) the ground state is degenerate 
independent of the parameters E and D. 

It is to be concluded that this degeneracy corresponds 
to a hidden symmetry of this class of Hamiltonian, akin 
to time reversal symmetry for the half-integer equivalent 
in zero field. This pseudo-time-reversal symmetry then 
implies a pseudo-Kramer's degeneracy. 



uncertainties and this minimum uncertainty state is un- 
squeezed. For finite E, 



h 2 S 
2 



D-E 



( h 2 S 



D + E 



D + E' y V 2 ) V D-E' 
Maximum squeezing occurs as E — > D. 



VI. QUANTUM FLUCTUATIONS - SQUEEZED 
STATES 



If S are not too small then S remains always close to 
the classical energy minima at = ir/2 and <fi = or 
7r, and to an excellent approximation near the first of 
these classical ground states Eqn. (||), after dropping a 
constant E° = -h 2 S 2 (D - E) and with h t = hg = 0, 
simplifies to, 



H=(D + E)p^ 2 + hpt + h 2 (D - E)S 



2/2 



(6) 



The field can be absorbed into the kinetic energy so 

H = (D + E) ( P<p + S) 2 +H 2 (D- E)S 2 2 

-(D + E)S 2 , (7) 

where S = h/2(D + E). Recall, S corresponds to displace- 
ment of S in m-space, however it is not equal to a s /2. In 
<j) space it adds a phase factor e zS( ^ to the wave function. 
Even for finite ft, this is just a harmonic oscillator with 
a ground state energy 



E u 



4^T^) +(n + ^ ' 



(8) 



where ujq = Sh\/D 2 — E 2 and where only the last term 
is of quantum origin. The wave function 



a{<t>) 



1/2 



D - E\ I 4XA q fD-^\l/2^ 

e l5(t) e- s y^+^) ^ . 9 
D+E ' v } 



Directly, the magnetic moment, 



(s z ) = (j>4>) = ns = h 



2{D + E)' 



(10) 



The same result follows from differentiating e n (ft). 

This ground state wave function of a harmonic os- 
cillator necessarily corresponds to the minimum of the 
uncertainty relationship AcfrAp^ > ft/2. This permits 
a natural definition of squeezed spin states. The trans- 
verse spin fluctuations are a minimum, e.g., for E = 0, 
a x =a v = h{S/2f/ 2 (e.g., a x 2 = (S x 2 ) - (S x ) 2 ). The S 
can be pictured as having fluctuating transverse compo- 
nents ~ ft(5'/2) 1 / 2 such that the total magnitude square 
is correctly h 2 [S 2 + 2 x (5/2)] = h 2 S(S + 1). Clearly 
the two quadrature components <j x and a y have equal 



VII. FORMULATION WITH THE \S Z ) BASIS 

The zero point energy associated with the transverse 
fluctuations is uoq/2 and correctly the low lying excited 
states have energies (n + (1/2))Hljq ~ h 2 SD above the 
ground state. However from Eqn. (|4|) the energy barrier 
is also ~ S 2 h 2 D and this would apparently correspond 
to a problem with a tunneling amplitude of order e~ s l 2 
which is not inconsistent with the Fourier transform of 
Eqn. (|j). However the unsqueezed states when E = 
are in fact eigenstates of S x , i.e., in this case tunneling is 
absent. 

It is necessary to carefully formulate the problem in 
terms of the basis set \m), i.e., in m-space. If \ijj) = 
Y^=-s a ( m )\ m )-> Schrodinger's equation in this basis 
can be written as (hi = h t =0): 

(e — (2Eh 2 m 2 — ftftm)) a(m) 
= \h 2 {D - E) \M^M2Xla{m + 2) 






M r ' 



M, 



(a{m - 2) 

Tl-l 2 ' 



}a(m) 



(11) 



,1/2 



where M\ = [S(S + 1) - st] l/z . The structure of this 
equation is worth noting. It has the form of two tight 
binding models, one for even m values and the other for 
odd m. The two chains are not connected as illustrated 



17/2 



13/2 




FIG. 8. Schrodinger's equation Eqn. (O) reflects two dis- 
connected chains. For whole integer S one chain corresponds 
to even and the other odd sites and are inequivalent. In con- 
strast for half-integer S the chains have n — ► — n symme- 
try. That these chains have the same eigenenergies reflects 
Kramer's theorem. 



in Fig. (8). This is not a periodic solid since (i) the 
diagonal energy increases steady going away from m = 
(ii) the chains have ends at m = ±S. This is rather a 
pair of discrete harmonic oscillator problems for large S 
when (ii) is not important since the wave function is well 
localized near m = 0. 

Making a continuum approximation, e.g., a(m ± 2) « 
a(m)±2(da/dn) + 2(d 2 a/dn 2 ), valid if S 2 (D-E) > 2E, 
and assuming large spin S ^> 1, this reduces to 



2Eh 2 - 2 



m 



e)a(ra) 



4h 2 S 2 (D-E) 



d 2 a(m) 
dm 2 



0. (12) 



This 



describes 
function 



a simple 

a(m) = 



harmonic 

-m 2 Q/2S 



oscillator with 
where, Q = 



\/(D + E)/(D — E) characterizes the degree of squeez- 
ing. This should be recognized as the Fourier transform 
of Eqn. (g). There is evidently no tunneling barrier and 
no tunnel ground doublet splitting! However there is a 
degeneracy since this same equation results for both the 
odd and even site chains. The splitting of these two states 
disappeared upon making the continuum approximation. 
As will become evident below, the discrete problem de- 
scribed by Eqn. ( |ll| ) has asymptotic corrections which 
differ as the m arguments are even of odd. It is an inter- 
esting observation that the discrete harmonic oscillator 
and tunneling problems are equivalent. 

The wave function for the eigenstate \S X = S) is 

-.1/2 

a°(m) = I 
lution to Eqn. 



(s- 



(25)! 



1/2 



and a°(m) has to be a so- 
E = 0, this corresponding to 



-m)\(S+m)\ 

Ll|) with h 

the absence of tunneling. Using Stirling's approximation 
this reduces to a(m) oc e~ m Q/ 2S with E = 0, however 
the deviations from this approximation are very impor- 
tant. The absence of tunneling when h = E = is put in 
evidence by a transformation: a(m) =a°(m)f(m). After 
some algebra Schrodiner's equation for f(m) reduces to: 

(e - (2Eh 2 m 2 - hhm + E )) f(m) = \h 2 (D - E) 
x '(S(S - 1) + m 2 )(f(m + 2) + f(m - 2) - 2/(m)) 



•n(2S-l)(/(m + 2)-/(m-2)) 



(13) 



When h = E = this evidently admits the solutions 
f(m) = 1 and f(m) = e l7Trn corresponding to the degen- 
erate states \S X = ±S). 



VIII. INTERMEDIATE SPIN FORMULATION 

With again the definition f((j)) = J2m e* m( ^/(m) this 
yields the exact Schrodinger's equation: 

(£ - 2Ep^)f(4>) = V(4>, P<p )f(<P) (14) 

with 



V{4> )P<t> ) = \{D-E)x 
2 sin 2 <j>(H 2 S(S -l)+p<p 2 )- ih(2S - 1) sin 2#^, 



(15) 



where £ = (e — E°). This is of the intermediate spin form 
but with a momentum dependent potential, V(0,p^). 

The effect of the transverse field h is twofold (i) near 
the classical ground state the spin sees a longitudinal 
field S(D — E) along with the transverse h and this tilts 
the axis of quantization towards the z-direction. (ii) The 
field h modifies the tunneling amplitude. 

In accounting for (i) the term 2Ep^ 2 is not relevant 
and so what is required is a solution of, 



(£ + h P4> )f(<J>) = V(<J>,p <p )f(<J>). 



(16) 



Needed is a small rotation about ?/-axis. The new wave 
function is defined by f((j)) = (l — (a s p ( f ) /2hS))g((j)) where 
h = h(D — E)cx s /2 1 this admixes, e.g., a little of the first 
excited into ground state. For large 5, g(4>) obeys 



_a s 



£g(<l>)=V(<l>,p 4 , + h-±)g((l>), 



(17) 



i.e., the field has been replaced by the displacement ^f. 
In order to account for (ii), a similar displacement is 
included in the full equation. If g{(f)) is the solution of 



(S - 2E ( P4> + h^f + h 2 E < ^-)g{ ( j ) ) 



a s 



= V(<J>, P<p + h^-)g(<J>), (18) 

then f((j)) = (1 — (a s p ( f ) /2hS))g((f)) is the solution to 

(£ - 2E P<p 2 + 2h(D + E) ^) f{4>) 

= V(4>,P4>)f(4>). (19) 

Following the discussion of intermediate spin, the shift 
determines the wave vector k = a s /2. In turn the 
statistical parameter must be chosen to give ground 
state energy determined without tunneling, i.e., Eqn. (g). 
Required is that h 2 E^f- = 4( ^ +£ n , so k = ^ = 



h/2Hy^2E(D + E). This fixes the Zeeman term in Eqn. 
©; 2H(D + E)^p^ = y/(D + E)l2Ehp4>. Since, 
in general, this is larger than hp^ the solution in the 
non-singular gauge A(l — (a s p ( f ) /2hS))uk((j)) + B(l + 
(a s p ( f ) /2hS))u-k((j)) involves a linear combination of both 
bare fields ±<\/(D + E)/2Eh such that the expectation 
value is equal the physical applied field. This requires 
A = 1 + ^2E/(D + E) and B = 1 - ^2E/(D + E). 
The total energy for level n and wave vector k is, 



E n (k)=E° 



h 2 



4(2? + E) 



+ [n+-) hiuo + e n (k). (20) 



The first two terms are purely classical in origin while 
the third corresponds to the harmonic motion about the 



bottom of the well. It is the last term which reflects the 
tunneling between wells. 

The problem therefore reduces to the determination of 
e n (k) and when the tunneling approximation is well jus- 
tified this reflects the tight binding approximation. Ac- 
tually for the case consider so far, with only a transverse 
field ft, the problem has a higher symmetry, see Fig. (9). 
The potential V((j),p<i>) has a <p period of a = ir rather 
than 27r so the reciprocal space unit vector K = 2. Since 
K corresponds to the smallest relevant displacement in 
m-space, it follows that a given solution only involves m 
values which differ by 2, as already noted in Sec. VII. 
The tunneling energy, 



of m are also half-integer and thus the solutions corre- 
spond to/c=± + ^f = ± + (h/2h^2E(D + Ej) and 
k = -^ + if = -\ + (h/2h^2E(D + E)) which yields 
a splitting, 



A = Aq sin ■ 



nil 



(23) 



£n(k) = COS7r/c 



(21) 



where Ao/2 is the yet undetermined (diagonal) matrix 
element for tunneling for level n. The evaluation of this 
will be addressed in Sec. XIII. There are two solutions 



2h^2E(D + E) 
i.e., the cosine is replaced by a sine. 

X. CONTINUATION OF SOLUTION 

The solution described above is based on a small ft ap- 
proximation it is valid in the very large S limit. That the 
matrix elements of h p — ph\\ fti = 2h^/2E(D + E) to 
excited states be small requires 2Ep < S 2 (D — E). For 
the first non-trivial whole integer point when ft = h\ 
this is equivalent to the requirement that the contin- 



which are acceptable for ft = 0, namely k — %£■ = uum approximation be valid, or that the shifts ±^f, or 
corresponding to even m values and k = ¥f- = 1 and 
odd m. Clearly a shift by one brings the set of integers 
m-values back to itself. These two k values correspond 



to e(k) = ±^f* and so the zero field tunnel splitting is 
reflected by the parameter Aq . For finite transverse fields 
ft the corresponding k = %*■ = (h/2hy/2E(D + E)) and 
k = 1 + ^ = 1 + (h/2Hy/2E(D + E)) and so the tunnel 
splitting is, 



A = Aq cos 



7rft 



2hyj2E(D + E) 



equivalent to a result first obtained by Garg 
tributed to topological quenching. 



(22) 



and at- 



IX. HALF-INTEGER SPIN 

Nowhere in the development has use been made of 
the assumption that S is an whole integer. The only 
difference when S is a half-integer is that the values 



^V(<P) 




^<p 



vels 



FIG. 9. When ht = he = the symmetry is higher and the 
period is tv rather than 2tt. The states with label n = 0, 1, . . . 
near the bottom of the wells are to a good approximation 
those of a harmonic oscillator. The box contains the unit cell, 
however in the singular gauge it is useful to continue the wave 
function and potential into the entire interval {— oo, +oo}. 



±fc, for a field ft = h p are smaller than the width of 
Gaussian wave function e~ m Q/ 2S . The wave function 
A(l - (a s p /2ftS)H(0) + B(l + (a s p^/2hS))u- k ((l)) is 
therefore a valid approximation for ft = fti for not too 
small S. However it will fail for modest values of 2E 
and large fields ft. The theory can be extended to cover 
such larger fields by a process similar to analytic contin- 
uation. The current approximation can always be devel- 
oped to construct the essentially exact solution a 1 (?n) for 
ft = fti. The dimensionless moment mi ~ h/2h(D + E) 
for fti can be determined. Re-defined are the states 
\<t>) = Em e^ m - m ^\m} so that S z = -ih{d/d<j>) + ftrai. 
This translates back to the origin the wave function 
a}(m). Using the original a°(n) an equation for f{n) 
is developed. Important is that in this equation there 
can be no term linear n reflecting an effective field, since 
this would imply a finite displacement ^ = k and hence 
a finite first derivative of the tunneling energy. This con- 
tradicts the assumption that the initial ft = fti reflects a 
whole integer point and hence a maximum in the tunnel 
splitting. The rather obvious point being made is that 
the whole (and half-) integer points can be represented 
without a mixture of different effective fields while in be- 
tween these points this is not possible. Continuing the 
development for points near to ft p , the new wave vec- 
tor k = 2f = 5h/2h^2E(D + E) where Sh = ft - h p . 
(The period cannot change since otherwise the current 
solution would predict a whole integer point which is not 
exactly at zero field.) Always the shifts are smaller than 
the Gaussian width and the new ^f differs from the old 
value by an integer and is hence equivalent. The fact 
that there is more than a single fashion by which to con- 
struct the wave function reflects the fact that each of the 
sets |0) independent of mi are (over-) complete. Impor- 
tant is the observation that away from the whole and 
half-integer points the wave function reflects an essential 



mixture of fields. The description with a minimum of dif- 
ference has fields h p ±^(D + E)/2ESh; 5h = h — h pi and 
where h p corresponds to the closest whole or half-integer 
point. 



XI. EFFECT OF OTHER FIELD COMPONENTS 

Consider first the effect of hi alone. This adds a po- 
tential, 



V((t>, p<f>) => V(<j>, P4) + HShe cos ( 



(24) 



and removes the symmetry between cf) = and it. The 
period of the potential is now the full a = 27r, with K — 1, 
but with two wells per unit cell. The harmonic levels near 
(j) = 7r have quantum numbers designated by n and are 
higher than those near </> = and which will be labeled 
n' '. Ignoring tunneling these have energies, 

h 2 ( 1\ 

E n =E - 7rR ; —^ + hnh e + In + - fio; , (25) 



4(L> + £) 



and, 



£^r? 



En 



tun! hi 



- J &jo. (26) 



4(1? + £) 

For a given /i there is only one solution with k = %£- = 
(h/2h^/2E(D + E)), the solution displaced by unity be- 
ing the same state. There are now two possibilities (i) 
E n y^ E n f and there is a very narrow band formed about 




FIG. 10. The effective potential ht — 0, hg ^ has two in- 
equivalent minima per unit cell. While the wells near ±7r are 
physically equivalent, in the singular gauge the wave function 
is multi- valued and so the phase when <j) increases by 2tt cor- 
responds to the Berry phase 2nk and is relevant. When the 
field h — the wave functions can be made real. Illustrated 
are the two different resonant cases. Above is the case when 
n — and n' — 2, i.e., when both wave functions are even 
and there is no change in sign. In contrast, below n — and 
n = 1, i.e., one wave function is odd and so the wave must 
change sign as its <j> argument advances by 2tt. 



both these energies. The energy differences have a very 
small dependence on k and hence h, see below, (ii) When 
E n « E n i there is resonant tunneling. There are again 
two cases (iia) when both n and n' are odd or even and 
(iib) when one of n and n' is odd and the other even. 

As illustrated in Fig. (10) there are two tunneling ma- 
trix elements one A n ^ n ' for the barrier between = 
and 7r and A n / >n for that between = — n and (or tt 
and 2tt). With first h = (h t =)0, Schrodinger's equation 
admits real solutions, ifj(m) or i/j ((/)), and so these matrix 
elements can be made real. For case (iia) when the sym- 
metry of the wave functions is the same the tunneling 
matrix elements A n>n / = A n / ?n = Aq. For finite /i, for 
the Floquet (Bloch) wave function, the phase advances 
by 2i\k from one cell to the next while to agree with 
the smaller unit cell appropriate when hi = the phase 
should advance by i\k from one well to the next within the 



cell. This leads to A r 



Aq and A^ 



C An. 



At this point it is possible to introduce the second 
transverse field h t by slight of hand. It is observed that 
if E changes sign in the original version of 7Y, Eqn. (g), 
then the two transverse axes change role and thus h be- 
comes the equivalent of h t . If T~L is considered to be a 
function of a complex parameter E the solution might 
be analytically continued from positive to negative val- 
ues. In particular k = (h/2h^/2E(D + E)) becomes 
pure imaginary, i.e., k - 
and it follows that A n?n 

e -27r(fc+zfct)^ o This information can be represented by a 
two level model: 



ih = i{h t /2hyj2E(D - E)) 

= e i<k+ik t ) AQ with An/ ^ = 



H = (E n 



'E n r)S z + — 



T(k+ik t ) C+ I e —i7r(k+ik t )Q- 



where the center of gravity energy \{E n + E n >) has been 
dropped. Degenerate perturbation theory then yields, 



£k 



±2V(^n 



at resonance 
h 



and kt 



E n ,y + 4A 2 | cos^(fc + ik t )\ 2 (27) 

±-^1 cos7r(/c + ikt)\, with k — 





The hi = h t 



2hs/2E(D+E) ~ 2h^/2E(D-E)' 

result is recovered when n = n' ' . That this is a good 
approximation is illustrated by comparison with exact 
results in Fig. (4). When for (iib) one wave function 
changes sign while the other does not it is necessarily the 



case, for h — h t = 0, that A n > 



-A n / 



Aq whence, 



Cfc 



±\\l(E n 



E n ,) 2 + 4A 2 | sin7r(/c + ik t )\ 2 (28) 



and 6/c = ±-^1 sinir(k+ikt)\ at resonance and, Eqn. (|22 
this is equivalent an original half-integer S. For case (i), 
Eqns. ( p7pg| ) serve to illustrate that when \E n — E n >\ > 
2Aosin7r/c, e.g., the upper sign of Eqn. (27) gives e^ ~ 
(E n — E n t) + 2(Ao 2 sin 2 irk/(E n — E n i)) which contains 
a /c-dependent correction which is smaller than those for 
cases (ii) by a factor of Ao/(E n — E n >). This factor is 
usually extremely small for the region between crossings. 



XII. COMBINED INTERMEDIATE SPIN 
EFFECTS 



determined for h = 0. The matrix element of 2ES X 
between the states m and m + 2 is 



Both h and hi cause the effective spin 5 - 



to alter- 



nate between effective whole and half-integer spin and 
the two fields can be combined so that, e.g., a field 
h = h^/2E(D + E) transmutes a whole into a half- 
integer spin. Adding hi = h\/D 2 — E 2 brings say the 
n = level in the higher well into resonance with the 
n = 1 of the lower well and the system again behaves as 
whole integer with a maximal tunnel splitting. However 
the fashion in which the two fields work is different. The 
transverse field h leads to a smoothly changing (Berry) 
phase shift of 5 = ir(hh/ \/2E(D + E) in the boundary 
condition ip((f) + 27r) = e z5 ^((/>) so that values of h which 
are not multiples of h^/2E(D + E) correspond to an in- 
termediate spin value. In contrast the phase shift caused 
by hi results from the changes in sign of the wave function 
for the odd n levels, i.e., when this field brings two levels 
into resonance 5 = 7r(n — n') which causes the effective 
spin value to alternate between whole and half-integer 
without a meaningful intermediate case. 



XIII. THE TUNNEL SPLITTING. 

There are three regimes for which different approxi- 
mate methods must be used in order to calculate the 
tunnel splitting, (i) If E <C D, and for quantization 
along the easy axis, the maximal tunnel splitting Ao is 
calculated treating E as a perturbation, see below, (ii) 
If S 2 (D - E) < 2E, i.e., very close to the point E = D 
the energy levels can be calculated using the hard axis 
for quantization and (D — E)S y is a weak the pertur- 
bation. To a very good approximation the energies are 
2Dh 2 m 2 - (h 2 /16E) where m = n+ (h/AHE) and n is 
an integer [||. Thus, e.g., the energy between the ground 
state and first excited state for h = is 2D, i.e., very 
large and nothing to do with tunneling. As a function 
of field the splitting is a sawtooth. Otherwise (iii) for 
E « D the problem corresponds to a dominant hard axis 
anisotropy with again (D — E)S y small. This is the case 
considered in considerable detail by Preda and Barnes 
[ |l2| . Following a method described by S. Coleman O] 
it is possible to calculate the tunneling matrix elements 
near to any crossing ra, ra' '. The result is |l2[ , 



Ao 
4 



= hujQ- 



^0 — + 
a — L rt 



n+1 



V K \\ 



) " exp (_2V2^p) . (29) 



For the remaining case (i) the energy levels are E m = 
-Dh 2 m 2 ] m = ±5, ±(S - 1), ... ± 1, 0. There are two 
cases (a) when the two levels in resonance n and n' both 
have the same parity and (b) when the parity is different. 
The whole integer case (a) is simplest since Aq can be 



n ~ — h 2 M rn + 1 M m + 2 



(30) 



recalling that M^ +1 = [5(5 + 1) - m(m + 1)] 1/2 . Thus 
the expression for the tunnel splitting is, 



. =(^y 



Mr +1 Mri 2 - 



M n+3 M n+4 

M n 'r^M n 'r 2 

lv± n' — A ±VJ -7i' - 



L n'—4: lv - L n , — 3 E ,_ 



M»,Z2M", 



ri-A ±v± ri-2> E n ,_ 2 -E 1 



-K>-X,_ v 



(31) 



For a half-integer point, i.e., case (b) it is observed that, 
with k = h/2h^2E(D + E), Eqn. ^j gives a splitting, 



5(h) = ttA 



2h^2E(D + E) ' 



(32) 



for small enough h. This must agree with a perturbation 
theory in which h can be taken to the smallest quantity. 
With the current axis of quantization the term hS x im- 
plies a matrix element of hM™ +1 /2 between the levels m 
and m + 1. The field provides one such matrix element 
to the expression, 



m=n-\-l 

M n+l M n+2 

. . . ±v± n ±v± n _^^ 



V 2 J 



M n+3 M n+4 

-^n+2 — -tin At,+4 — ^n 



171—1 771—1 



1 



M m+1 



E n 



1 



E, 



'm+l 



E r , 



M m+2 M n+3 

m m+l m n+2 ' ' • 



M n ~ 3 M n ~ 2 

ti n ' -2 — ti n 



' Miz\Mi,_ x 



= C(n,n')h 



2(D + E) 



(33) 



so that, 



A = C{n,n') — 



E 



2(D + E) 



E 



C(n,n') 



E 



2{D + E) 



(34) 



which has the same behavior as Eqn. ( pl|) for transitions 
between the same parity. In both cases, 



A ^2tl 2 Ee- i;i ^ 2LlSo ; S = \n 



E 



2(D + E) 



(35) 



i.e., 2E is the attempt frequency and the action agrees 
with that which is obtained for the ground state to 
ground state tunneling using the functional integral 
method |J . There is little interest in reducing the expres- 
sions Eqns. (31) and (|3J) since precise values for splitting 
are easily obtained numerically. It is an interesting result 
that the attempt frequency is 2h 2 E and not perhaps the 
characteristic energy 2h 2 ^2E(D + E) or h 2 \/D 2 - E 2 . 



XIV. OTHER SYMMETRIES 

The issue that intermediate spin effects are manifested 
even when there is no evident easy plane has been raised 
in the Introduction. Consider cubic symmetry, the ap- 
propriate Hamiltonian, 



H = a [S x 2 + S y 



hS z 



(36) 



has eight classical equilibrium directions, Fig. (lc) which 
with a positive do not coincide with the axes. Quan- 
tum mechanically, in the absence of a field the system 
resonates between each of these equilibrium directions. 
However in the language of Sec. VII, there are only four 
chains and for large spin only four closely spaced levels. 
Symmetry dictates that in zero field the chain with sites 
...,—3, 1,5,... has levels which are degenerate with the 
similar chain with n — > — n. To a good approximation 
for large S Schrodinger's equation is: 



[e — (ah 4 m 4 



3afr 4 



(S(S+l)-m 2 ))]a(m) 



TT E(S(S + 1) - (m ± 2) 2 )a(m ± 4) 



(37) 



where a(m) is the wave function. An analytic approach 
to the solution of this equation is facilitated by a modest 
field h and large S. This suppresses the tunneling be- 
tween the state with positive and negative projections on 
the z-direction. Now tunneling occurs between the four 
classical equilibrium directions with a negative S z compo- 
nent. It is the four least action paths between these which 
map to a circle. Making the continuum approximation 
near the minimum in the the diagonal energy with a neg- 
ative such projection defines a discrete simple harmonic 
oscillator with a potential ~ 2ah A S 2 (m — mo) 2 where 
mo is the equilibrium dimensionless magnetization. In 
Fourier space this becomes the kinetic energy and the 
potential ~ afi 4 5' 4 cos4</> so that there are four equiva- 
lent minima. Apart from this reduction of the period 
by a factor of two, the problem is then fully equivalent 
to that considered in Sec. VIII and the tunnel splitting 
quasi-periodic as illustrated in Fig. (6). As the field in- 
creases the minimum energy directions bend towards the 
z-axis and this reduces the potential barrier between the 
equivalent wells and increases the kinetic energy, hence 
the steady increase in the tunnel splitting. The separa- 
tion between the half-integer points when the splitting is 



zero is seen to be ~ Ah A S A a but varies since the coeffi- 
cient of the m 2 kinetic energy term also increases with 
field. In fact there are two level crossings at small fields 
which occur before the large field approximation is valid 
and which are not visible in this figure. At least for larger 
fields these half-integer points really do correspond to a 
Kramer's degenerate half-integer spectrum for the low ly- 
ing quartet as is explained in the Introduction and the 
Fig. (6) caption. 

The symmetry appropriate to Mni2 has an easy z-axis 
with equivalent x and y directions so that the Hamilto- 



nian 



H = -DS Z 2 + a [S x 2 + S y 2 ] + hS x , 



(38) 



where, with a positive, the x-directed transverse field is 
along a hard direction in the x — ^/-plane. Classically 
in zero field there are equivalent equilibrium directions 
which have the magnetization parallel or anti-parallel to 
the ^-direction. A finite h does no alter the equivalence 
of these two directions, the magnetization simply gains 
a component opposite to the field direction. However 
the field does have an effect on the tunneling paths. In 
zero field there are four least action paths which pass 
via one of the easy directions in the x — ?/-plane. With 
a field in this particular direction and a positive there 
remain only two equivalent paths which pass via an easy 
direction for which the magnetization has a component 
in the direction away from the field. It is these two least 
action paths which map to the intermediate spin circle. 
Here the paths are not co-planar and the field h will enter 
in the potential but without altering the symmetry. 

It is possible to analytically account for all four tunnel- 
ing paths for large S and not too large a h by imagining 
that there are two crossed circles. A given circle com- 
prises two easy directions in the x — y-piane and z-axis, 
i.e., the two planes containing the circle make an angle 
of 45° to the applied field. The wave function for each 
circle is calculated and the two are added to give the fi- 
nal result. Because the field is at 45° to each circle the 
wave vector has equal real and imaginary parts and the 
net result is that displayed in the caption to Fig. (5). 



XV. EXPERIMENTAL VERIFICATION OF 
FIELD MIXING 

Certainly the most striking claim is that the wave func- 
tion A(l-(a sP(f) /2hS))u k ((l)) +5(1 + (a a ptl2hS))u- k {<l>) 
for an intermediate spin system is a mixture of a 
state with the local applied field of magnitude h p ± 
sf{D + E)/2ESh; Sh = h — h p , i.e., fields both larger 
and smaller than the distant applied field h. (With an 
equivalent claim for ht which involves yiD — E)/2Ehf.) 



10 




FIG. 11. The dimensionless moment m(h). The field is 
divided by half the period so that all steps coincide and S = 5. 
The moment period is exactly unity for E — D but otherwise 
is less and implies a Schrodinger's cat with a superposition of 
fields. 

For a particle on a ring this would imply a distortion in 
the near field with a range of order of the ring diam- 
eter. Unfortunately the spin has no natural dimension 
and the ring diameter must be scaled to zero. Exper- 
imentally the existence of this superposition is easy to 
establish. In the absence of the superposition the wave 
function would have a simple Floquet form: e lk< ^ Uk{(j)) •> 
where k = ^f is determined, modulo unity, from the ex- 
perimental tunnel splitting: Ao|cos7rfc|. Theory gives 
k — %£- = Sh/2h^2E(D + E) which agrees with exper- 
iment for Feg. Following the Fourier transform to 
m-space the factor e lk< ^ corresponds to a displacement of 
k and would lead to a dimensionless magnetic moment, 



^non-mix 



k ~ 2 



5h 



2hy/2E(D + E) ' 



(39) 



again modulo unity. In the limit of large S the theory 
accounting for the superposition leads to 



m n 



2h(D + E) 



D. This has yet to be 



which is smaller except when E 
verified experimentally. 

A more subtle way of stating the same principle goes 
as follows: The transverse field h couples to the physi- 
cal S z , i.e., this operator in the non-singular gauge, and 
which has a spectrum with a period h or in units of mag- 
netic moment g\iB^ where g is the g- factor. Thus if M z is 
the measured moment in the direction of the transverse 
field, m = M z /giiBh, the dimensionless magnetization, 
has a natural period of unity, a statement equivalent to 
the observation that physical effects are periodic in the 
flux quantum <I>o- This would be the case if the state did 
not consist precisely of a superposition of field directions. 
Because of the admixture, the measurable period is re- 
duced to Am = mi = J p^ E < 1 as illustrated with ex- 
act results in Fig. (11). It should be noted that the value 



of mi given above is not very accurate almost uniquely 
because dimensionless magnetization m > h/2h(D + E) 
for any reasonable S because of admixture into excited 
states. For example with E/D = 0.35 and S — 10 it 
was found that m = 1. 0353718 {h/2h{D -\- E)) was indis- 
tinguishable from the numerical data for m in the range 
< h < hi and gave with five figure precision the numer- 
ical mi. (Notice m\Kh = 4E while (hh/rrti) = 2(D+E).) 
It is also possible to argue for the local field 
y/{D + E)/2Eh in terms of the two level model. For an 
integer spin system for the ground state with h = hg = 
this is TL = 2ShiS z + AqS x where, as above, Ao is the 
tunnel splitting. For half-integer spin, following Kramer's 
theorem this slitting is absent. However since the trans- 
verse field breaks time reversal symmetry it should create 
a splitting. The only satisfactory fashion of introducing 
this transverse Zeeman effect without introducing a new 
parameter is to write TL = 2ShgS z + A§mS x , where m 
is the dimensionless magnetic moment in the x-direction. 
Since m = h/2h(D+E) this gives a splitting 2(A /H(D+ 
E))h. However from Eqn. ( p3[) this is correct only if h 
is replaced by the local applied field tt^/(D - 



to give a splitting irA Q (h/2h^/2E(D + E)). 



f E)/2Eh 
The field 



h t similarly generates a term 7rA(h t /hy/2E(D — E))S y . 
Clearly this model it is subject to experimental verifica- 
tion and nicely illustrates the effects of squeezing via the 
different factors involved with the different components 
of the transverse field. 



XVI. EXPERIMENT AND SPIN SQUEEZING 

The lack of symmetry between the x and ^/-directions 
reflects spin squeezing. For equal transverse fields 
m x /m y = y/(D-E)/(D + E) = cr x 2 /a y 2 = Q, the de- 
gree of squeezing. More directly, e.g., a 2 is given by sum- 
ming the intensity of all transitions measured by the lon- 
gitudinal response to an r.f. field in the x-direction. As 
explained above, the transverse matrix elements within 
the ground doublet are ~ Aom, i.e., extremely small. It 
follows that the transitions with energy ~ D exhausts 
this sum rule. Experimentally neutron scattering which 
can measure the same correlation function in zero static 
field would seem most appropriate. 



APPENDIX 

Given that h 2 D is the dominant anisotropy energy it 
would seem the most logical to use this easy axis for the 
quantization, however much of the physics of the prob- 
lem corresponds to the fact that term proportional to E 
creates an easy plane. It is therefore useful to re- write 
TL quantized along the hard axis. To this end, first use 
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S y 2 = h 2 S(S + 1) - S x 2 - S 2 whence, 

H = -(D - E)S 2 + 2ES 2 - Eh 2 S(S + 1), 

then interchanging the x and z axes results in 

H = 2ES 2 - hS z -{D- E)S X 2 - h £ S x - h t S y , 

where the constant has been omitted and the fields re- 
introduced. 

The case when D < E can be transformed to D' > E' 
where D' = 3E - D and E' = D + E. To show this, first 
use again S z = H 2 S(S + 1) — S x — S y to write: 

H = -DS(S + 1) + (D + E)S X 2 -(E- D)S y 2 . 

With the interchange of the z and ?/-axes and to within 
a constant: 

H = -(E- D)S 2 + (D + ^)^ 2 . 

Comparing this with W = -(25 - ^)^ 2 + 2ES 2 it is 
seen that, 
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D + E^»2E, E-D^D-E 

reduces the second to the first version. These manipu- 
lations leave the term hS z unchanged. Hence it must 
be, in particular, that the result for the period h\ = 
2Hy/2E(D + E) be invariant under the transformation 

E ■ 



\{D + E) D^ l -(ZE-D) 



Then if E > D it follows that D f = 3E - D > 2E 
while E' = D + E < 2E so that for the transformed 
problem D 1 > E' , i.e., only that case when D > E need 
be considered. 

Another useful version of TL is for h t = hg — 0, 



n 



-(D + E)S X 2 -2ES 2 -hS z . 
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